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ABSTRACT 

I propose a mechanism for axisymmetrical mass loss on the asymptotic giant branch (AGB), 
that may account for the axially symmetric structure of elliptical planetary nebulae. The proposed 
model operates for slowly rotating AGB stars, having angular velocities in the range of 10“^WKep ^ 
10 ^WKep, where WKep is the equatorial Keplerian angular velocity. Such angular velocities 
could be gained from a planet companion of mass ^ O.lMjupiter, which deposits its orbital angular 
momentum to the envelope at late stages, or even from single stars which are fast rotators on the 
main sequence. The model assumes that dynamo magnetic activity results in the formation of cool 
spots, above which dust forms much easily. The enhanced magnetic activity toward the equator 
results in a higher dust formation rate there, and hence higher mass loss rate. As the star ascends 
the AGB, both the mass loss rate and magnetic activity increase rapidly, and hence the mass loss 
becomes more asymmetrical, with higher mass loss rate closer to the equatorial plane. 

Key words: Planetary nebulae:general — MHD — stars: AGB and post-AGB — stars: mass loss — 
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1. INTRODUCTION 

The structures of planetary nebulae (PNs), among other arguments, suggest that many AGB 
stars, the progenitors of PNs, have a high mass loss rate at the end of their AGB phase, in what 
is termed superwind. The mechanism responsible for the high mass loss rate is thought to be 
radiation pressure on dust (e.g., Knapp & Morris 1985), with some contribution from molecule 
acceleration (Johnson, Alexander, & Bowen 1995). The rapid growth of the mass loss rate at the 
end of the AGB has been suggested to result from the increase of the density scale height above 
the photosphere (Bedijn 1988; Bowen & Wilson 1991). Since the density scale height H appears 
in the argument of an exponential term in the expression for the density, p{r) oc both the 

density at the location of dust formation and mass loss rate increase rapidly as H increases. 

In many elliptical PNs the inner shell, which was formed from the superwind, deviates more 
from sphericity than the outer region, which was formed from the regular slow wind (prior to the 
onset of the superwind). In extreme cases the inner region is elliptical while the outer region (outer 
shell or halo) is spherical (e.g., NGC 6826). In addition, most (~ 75%) of the 18 spherical PNs 
(listed in Soker 1997 table 2) do not have superwind, but just an extended spherical halo. The 
correlation between the onset of the superwind and the increase in the asymmetry of the wind is not 
perfect, and in some cases the inner and outer region both have a similar degree of asymmetry (e.g., 
NGC 7662). The correlation between the superwind and highly non-spherical mass loss, although 
not perfect, is a challenge to any model for the formation of elliptical PNs and for the superwind’s 
high mass loss rate. 

In an earlier paper (Soker 1995) I proposed that this correlation results from a late interaction of 
the progenitor AGB star with a binary companion, in most cases planet or brown dwarf companions 
(Soker 1997). As the star expands along the AGB the ratio of its radius i?* to the orbital separation 
a increases. Since the tidal spiraling-in time goes as (a/i?*)®, there is a rapid transition from 
negligible tidal interaction to strong tidal interaction. The tidal interaction results in the spinning 
up of the envelope, and in most cases the companion enters the envelope and forms a common 
envelope. The rotation of the envelope results in both axisymmetrical mass loss and a higher mass 
loss rate. This model requires that the companion be at the right orbital separation range; it 
cannot be too close in order to avoid interaction when the star is on the red giant branch (RGB), 
and it cannot be too wide, otherwise no tidal interaction will occur. Main sequence stars of mass 
Mms ^ 2Mq expand to large radii on the RGB, and therefore are likely to interact with their 
companion already on the RGB (Soker 1998). Therefore, the late tidal interaction model of Soker 
(1995) can be applied only to stars having main sequence mass of Mms ^ 2Mq. 

In the present paper I propose a different mechanism to explain the correlation of superwind 
and axisymmetrical mass loss. In this model I assume that a weak magnetic field forms cool stellar 
spots on the surface of the AGB star. The cooler spots facilitate the formation of dust closer to the 
stellar surface, enhancing the mass loss rate (Frank 1995). Enhanced dust formation above cool 
regions of red giants was assumed by Schwarzschild (1975), where in his model the cool regions 



- 3 - 


are formed by very large convective elements. Frank (1995) conducted a detail study of dust 
formation above cool AGB starspots, and showed that the mass loss rate above the spots increases, 
though the terminal wind velocity does not change much. Frank (1995), though, does not discuss 
the source of the cool starspots. A slow rotation, which is required for the amplification of the 
magnetic field, forms a dipolar field which has a stronger magnetic activity near the equatorial 
plane; hence the axisymmetrical mass loss. One advantage of the proposed model is that only a 
slow rotation on the AGB is required, which can be acquired from a gas giant planet, or even by 
singly evolving stars (i.e., no further spin-up during the post main sequence evolution is required), 
if they are rapid rotators on the main sequence. The model proposed in the present study is built 
to explain the axisymmetrical mass loss from the progenitors of elliptical PNs. The model does not 
apply to bipolar PNs (PNs that have two lobes with an ‘equatorial’ waist between them), which 
seem to result from stellar binary systems (Livio, Salzman & Shaviv 1979; Morris 1987; Gorradi & 
Schwarz 1995; Soker 1998; Mastrodemos & Morris 1998; and more references in these papers). In 
§2 I summarize a number of previous papers dealing with dust formation and magnetic activity, 
which are important ingredients of the model. In §3 I present the phenomenological model, and in 
§4 I discuss the results. The summary is in §5. 


2. BASIC INGREDIENTS OF THE MODEL 

2.1. Dust 


Bedijn (1988) and Bowen & Wilson (1991) suggest that the rapid increase in the mass loss 
rate toward the end of the AGB phase results from the increase in the density scale height above 
the photosphere. The density scale height is given by oc R^T/M^, and the density profile is 
p{r) oc exp(—r/R), where is the stellar radius, M* its mass and T the temperature above the 
photosphere. The density at the place where the winds starts, in the case of radiation pressure on 
dust this being the dust condensation radius, is given by 

Pc - exp(-AM*/R*), (1) 

where A is a constant A few X 10^Mq^Rq (Bedijn 1988; Bowen 1988). The mass loss rate is 
proportional to this density. As a star ascends the AGB, R* increases and, because of mass loss, 
M* decreases, and both enhance H and hence the density pc and the mass loss rate. Bowen &: 
Wilson (1991) argue that dust formation is not required to explain the rapid increase in the mass 
loss rate, but it is required to obtain high mass loss rates as observed on the upper AGB. In the 
present study I assume that radiation pressure on dust is the mechanism responsible for the high 
mass loss rate. 
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Dorfi & Hofner (1996) proposed a model based on dust to explain axisymmetrical mass loss. 
Their model requires a large amount of angular momentum, despite what they term “slow rota¬ 
tion”. Although it is quite plausible that rotation together with radiation pressure on dust form 
axisymmetrical mass loss, I show now that their model must incorporate a binary companion to 
spin-up the envelope. In order for their proposed scenario to work, Dorfi & Hofner (1996) require 
the angular velocity of their AGB star, of radius R = 500Rq, to be ^ 10% of the Keplerian angular 
velocity. We hardly find main sequence stars with such high rotational velocity; there is no way 
a singly evolved AGB star can obtain this rotational velocity. Approximating the density profile 
on the AGB by p oc r^, where r is the radial distance from the center of the star, we find the 
envelope’s moment of inertia to be lenv = (2/9)Menv.R*, where Menv is the envelope’s mass. With 
this envelope density profile the angular momentum of the envelope decreases with mass loss as 
J oc (Harpaz &: Soker 1994). Since most of the mass loss occurs on the AGB, we find that the 
ratio of envelope angular velocity on the upper AGB to the Keplerian (critical) angular velocity for 
a single star evolution is given by 

/ ^ \ / -^MS \ ^ f ^env \ 

y O.lWKep/j^g \0.01i?AGB/ VO.lMenvO/ 

where the subscript MS means that the quantity is taken at the end of the main sequence, and 
MenvO is the envelope mass at the beginning of the AGB. The contraction of the star from the 
AGB during the post-AGB phase will result in a further spin down (Soker 1992). This is because 
the high mass loss rate causes mass shells to move outward in the convective envelope, despite the 
contraction of the star. Much more massive stars which retain massive envelopes may spun-up 
when they move to the left of the HR diagram during their blue loop (Heger & Langer 1998). This 
is caused by the transfer of angular momentum from the inner to the outer convective region, as 
the inner radiative region of the envelope expands (Heger & Langer 1998). Since the contraction of 
post AGB stars is due to the depletion of the very low-mass post AGB envelope by the wind, which 
carries most of the envelope’s angular momentum, this mechanism cannot work for contraction 
beyond the AGB (Soker 1992). 

The orbital separation of a low mass secondary when tidal interaction becomes significant is 
a ~ 5i?* (Soker 1998). If the secondary deposits all its orbital angular momentum to the envelope 
of mass Menv, the ratio of envelope angular velocity iv to the surface Keplerian angular velocity 
(jKep is given by 

~ 0.1 (—^) 

t^Kep VO.OlMenv/ 

assuming that the entire envelope rotates uniformly. We conclude that in order to spin-up the 
envelope as required by Dorfi & Hofner the secondary mass should be M 2 > O.OlMenv However, as 
Harpaz & Soker (1994) show, the envelope’s specific angular momentum of an AGB star decreases 
with mass loss as Lenv/Menv oc Mgnv Therefore, to supply the angular momentum for a longer time, 
the companion mass should be much larger than O.OIMq, i.e., a brown dwarf or a low main sequence 
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star. Other effects that such a companion can cause (Soker 1997) should then be considered as well. 
In the present paper the required angular velocity is much smaller, since I incorporate magnetic 
fields, and a planet companion of mass ~ O.lMjupiter is sufficient to spin-up the envelope. 

We now show the problems with direct effects of rotation, e.g., centrifugal force, for more 
general cases. The centrifugal force on the equatorial plane increases the scale height by a factor 
of [1 — (cu/wKep)^]”^- The scale height along the polar directions does not change. Inserting this 
factor in equation (1) for the equatorial plane we get {pc)e — exp — [ca/wxep]^)] • The 

ratio of equatorial to polar density is then 

~ exp[TM*i?7^(a;/wKep)^]- (4) 

[Pc)p 

The strong dependence of the density contrast, and hence mass loss rate contrast, between the 
equatorial and the polar directions on the angular velocity is what Dorfi & Hofner (1996) find in 
their numerical calculations. Numerically, the results of Dorfi & Hofner (1996), for M* = IMq and 
R^, = 500720, fit a value of H ~ 2 x 10'^Mq^Rq rather than A = 5 x 10 ^Mq^Rq found by Bedijn 
(1988). Two problems for a model based on the centrifugal force emerge from equation (4). First, 
as the star ascends the AGB and loses mass, M* decreases, 72* increases, and because of expansion 
and mass loss (cu/wKep)^ decreases. Therefore, from equation (4), the density contrast between the 
equatorial and polar directions decreases, contrary to what is observed in most PNs. The second 
problem is the required angular velocity. Inserting A = 2x 1O^Mq^ 720 in equation (4), we find the 
angular velocity required to obtain a density contrast q = {pc)e/{Pc)p to be 


OJ 


^Kep 


0.08 


/Ingy/^ 



72* 

200720 


1/2 


(5) 


As mentioned above, such a high angular velocity requires a companion to spin-up the envelope, 
at least a brown dwarf, and probably a low mass stellar companion. 


2.2. Magnetic Field 


Direct magnetic effects, through magnetic tension and/or pressure, have been suggested to 
determine the mass loss geometry from AGB stars (Pascoli 1997 and references therein), or to 
influence the circumstellar structure during the PN phase (Ghevalier &: Luo 1994; Garcia-Segura 
1997). As I show in this subsection, these models must incorporate a binary companion to sub¬ 
stantially spin-up the envelope. These models will be compared later to the non-direct magnetic 
effects proposed in the present study. 

In several papers, Pascoli (1997, and references therein) has suggested that magnetic activity 
plays a major role in determining the mass loss rate and geometry from AGB stars. Although 
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Pascoli’s general idea of magnetic activity might be right, there are a number of specific problems 
with the details presented in his models (Soker & Harpaz 1992). Pascoli (1997), for example, 
proposes that the magnetic field formation, through dynamo mechanism, occurs close to the core 
of AGB stars. He then assumes that the core surface, of radius Rc = 10® cm, rotates at a velocity 
of Vc = 100 km s“^, and that the envelope convection penetrates down to the core’s surface. Both 
assumptions seem to me problematic. First, detailed models of AGB stars show that there is a large 
radiative zone above the core (e.g., Soker 1992), and that the envelope’s convective region starts 
only at r ~ IRq. Second, the convection acts to enforce a constant angular velocity uj{r) ~ ujq 
throughout the convective region, contrary to Pascoli’s assumption of constant specific angular 
momentum (jj{r) oc r“^. Therefore, even if the convection penetrates down to the core’s surface, 
it will substantially reduce the angular velocity compared with the value assumed by Pascoli. 
Therefore, the rotation velocity assumed by Pascoli requires that the AGB star be spun-up by a 
companion. Further indication for the too high velocity assumed is the finding that single white 
dwarfs rotate very slowly, Urot ‘C 50 km s“^ (Heber, Napiwotzki, & Reid 1997). Using lower rotation 
velocities, as indicated above, will result in a much weaker magnetic activity than that assumed by 
Pascoli. 

The second model I find to require substantial envelope spin-up is that of Garcia-Segura (1997; 
see also Garcia-Segura et al. 1998), which is an extension of the model proposed by Chevalier &: 
Luo (1994; see also Chevalier 1995). This model is based on the tension of the toroidal component 
of the magnetic field in the wind: the wind in the transition from the AGB to the PN phase or the 
fast wind during the PN phase. Close to the star the magnetic pressure and tension are negligible 
compared with the ram pressure and thermal pressure of the wind. As the wind hits the outer 
PN shell, which is the remnant of the slow wind, it goes through a shock and slows down, and the 
toroidal component of the magnetic field increases substantially. This may result in the magnetic 
tension and pressure becoming the dominant forces near the equatorial plane. In particular, the 
magnetic tension pulls toward the center and reduces the effective pressure in the equatorial plane. 
According to this model (Chevalier & Luo 1994; Garcia-Segura 1997), then, the equatorial plane 
will be narrow, leading to an elliptical or bipolar PN. 

The efficiency of this model is determined by the parameter (Chevalier & Luo 1994) 

y M^Vw) ' J Ek V 'IJw J 

where Bg is magnetic field intensity on the stellar surface, rg the stellar radius, the mass loss 
rate into the wind, Vw the terminal wind velocity, and Urot the equatorial rotational velocity on the 
stellar surface. In obtaining the second equality the expressions for the magnetic energy luminosity 
Eb = 4:TrrgVwBg/8Tr and for the kinetic energy luminosity Ek = M^v‘^/2 were used. For the model 
to be effective it is required that a ^ 10“^, but a typical value of a ~ 0.01 is used by Garcia- 
Segura (1997). For the sun a ~ 0.01 and ~ 2 x 10“^ (Chevalier & Luo 1994). However, 

in the sun it is magnetic activity which determines the mass loss rate, as we see from the ratio 
Es/Ek — 500. It is commonly assumed that radiation pressure drives the winds of central stars of 
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PNs, and that pulsation together with radiation pressure drives the wind of AGB stars, and red 
giants in general. Therefore, the sun is not a good example of this model for singly evolved stars. 
If magnetic energy does not drive the wind then Eb/E^ ^ 1, and the model of magnetic shaping 
requires v^ot/vw ^ 0.01. Such rotation velocity is impossible for singly evolved AGB or post-AGB 
stars to attain (Harpaz & Soker 1994; §2.1 above). 

It is clear that in order for the magnetic shaping model of Chevalier & Luo (1994) and Garcia- 
Segura (1997) to be of any significance, a substantial spinning-up by a binary companion is required. 
But even if this condition is met, I find two other problems with this model concerning shaping on a 
large scale. The two problems result from MHD instabilities. Other MHD instabilities might exist 
as well (Livio 1995). (a) The magnetic shaping model requires that the magnetic field lines circle 
the central star in the equatorial plane. However, it is not clear that this will be the case. Because 
of MHD instability the magnetic field escapes from the sun in non-axisymmetrical magnetic flux 
loops (e.g., Bieber & Rust 1995; Caligari, Moreno-Insertis, & Schiissler 1995). Many flux loops 
which escape the sun do not circle the sun. It is possible, though, that when the magnetic pressure 
is low more flux loops will circle the star, (b) For the magnetic field to be of any significance, 
it should be regenerated by a stellar dynamo. The idealized toroidal magnetic field that results 
from a dynamo has opposite directions in the two stellar hemispheres (e.g., Bieber &: Rust 1995). 
Therefore, as the magnetic pressure becomes dominant after the wind slows down, I expect that 
reconnection of magnetic field lines close to the equatorial plane will occur. This will reduce the 
magnetic pressure and tension near the equatorial plane, and as a consequence will reduce the 
efficiency of the model. 

The last two problems relate only to the large-scale shaping. Magnetic field may still be 
strong but with a short coherence length. I think that when substantial spinning occurs and if 
dynamo activity becomes efficient, magnetic fields may play a substantial role in small-scale shaping, 
e.g., MHD instability modes on small scales, fragmantation of regions in which SiO masers exist 
(Hartquist & Dyson 1997), and non-thermal radio emission (Dgani & Soker 1998). 


3. THE MODEL 


3.1. Outline of the Model 


In order to circumvent the fast rotation required by the models described in the previous sub¬ 
sections, but still using the popular model of radiation pressure on dust particles for the superwind, 
I propose the following scenario. I assume a weak magnetic activity, i.e., the magnetic tension or 
pressure never become the dominant force globally, neither on the surface nor at large distances 
from the star. As in the sun, the magnetic activity results in the formation of cool spots. Above 
these spots, I suggest, dust forms much more easily, i.e., faster and closer to the stellar surface 



(Frank 1995). The formation of dust closer to the stellar surface means higher density there, and 
hence higher mass loss rate, as summarized in §2.1 and as was shown by Frank (1995). A rotation 
is required in this scenario to amplify the magnetic field, but a much slower rotation than what is 
required in the models mentioned earlier. It might even be the case that fast rotating singly evolved 
stars can possess such rotation velocities (i.e., no further spin-up is required, §3.2). Fastly rotating 
main sequence stars have mass of Mms ^ 2Mq. This is the mass range of progenitors of bipolar 
PNs (e.g., Corradi & Schwarz 1995), and of stars that form type I PNs (Torres-Peimbert &: Peim- 
bert 1997). However, not all massive stars form bipolar PNs, and not all type I PNs are bipolars 
(Corradi & Schwarz 1995). Soker (1998) estimates that only ~ 40% of stars having main sequence 
mass of 2Mq ^ Mms ^ 8Mq form bipolar PNs. If bipolar PNs are formed from interaction with 
stellar companions (Morris 1987; Livio 1997; Soker 1998), then most of the rest massive stars will 
form elliptical PNs. The massive progenitors of elliptical PNs may rotate fast when arriving to the 
AGB through interaction with planets (Soker 1996), or they may retain enough angular momentum 
from their main sequence phase. 

The dynamo amplification results in a higher activity closer to the equatorial plane, therefore 
a higher mass loss rate there. The increase in the mass loss rate contrast between the polar and 
equatorial directions as the star ascends the AGB is explained by the rapid decrease in the thermal 
pressure of the atmosphere, and the formation of dust at higher densities (§2.1). The cool spots 
mechanism means that the enhanced mass loss rate in the equatorial plane will be sporadic, though 
over a long time period it will be averaged to form a smooth nebula. Episodic, local formation of 
dust occurs close to the surface of R Coronae Borealis (RGB) stars (Glayton 1995 and references 
therein). RGBs are hydrogen-deficient carbon-rich supergiants, thought to result from the final 
helium shell flash (Renzini 1990). The common interpretation of observations is that dust forms 
close to the stellar surface (< 2R*), on an area covering only a fraction of the surface, and that 
it occurs at irregular intervals. RGBs are not rapid rotators and there is no evidence for binarity 
(Glayton 1995). Whitelock et al (1997) find erratic behavior similar to that of RGB stars in three 
Garbon-rich variable AGB stars. The ejection of carbon-rich dust occurs in a preferred direction, 
which they speculate is the equatorial plane. Frank (1995), whose model is built for carbon rich 
stars, brings further evidences to support clumpy mass loss from AGB stars. According to the 
model proposed here, magnetic activity, which requires only slow rotation, forms cool spots on 
which dust is formed. The formation of cool spots by magnetic activity is seen on the surface of 
the sun. The assumption of enhanced dust formation above cool regions of red giants was used by 
Schwarzschild (1975). In Schwarzschild’s model the cool regions are formed by very large convective 
elements. Another hydrogen deficient object is the PN A30 (Jacoby & Ford 1983). This PN 
has a large, almost spherical, halo, with optically bright, hydrogen-deficient, blobs in the inner 
region. The blobs are thought to result from a late helium shell flash. These blobs are arranged 
in a more or less axisymmetrical shape. It is possible that during the formation of the halo dust 
was forming far from the stellar surface. If after the helium flash the formation of dust occurred 
closer to the stellar surface, the process became more vulnerable to magnetic activity, resulting in 
the axisymmetrical mass loss. 
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A hint of the presence of magnetic fields in the atmospheres of some AGB stars comes from the 
detection of X-ray emission from a few M giants (Hiinsch et al. 1998). An even stronger motivation 
for the proposed model is the observation of a magnetic field in the extended atmosphere of the 
Mira variable TX Cam (Kemball & Diamond 1997). Kemball & Diamond (1997) find the intensity 
of the magnetic field in the locations of SiO maser emission to be ^ 5G. These local emission 
regions form a ring around TX Cam at a radius of 4.8 AU ~ 2R*. Kemball & Diamond (1997) 
mention the possibility that the mass loss occurs in a preferred plane. They also suggest that “The 
fine-scale features [of the Maser image] are consistent with local outflows, flares or prominences, 
perhaps coincident with regions in which localized mass loss has taken place.” 


3.2. The Magnetic Activity 


Because of the lack of basic theories for (a) dynamo operation in ACB stars, where the convec¬ 
tive overturn time is shorter than the rotation period, (b) formation of stellar cool spots through 
the action of magnetic flux tubes, and (c) dust formation close to the stellar surface, the required 
magnetic activity will be studied phenomenologically. For this I take the following: 

(1) From the sun it is known that the magnetic fields in sunspots (Bgpot ~ 2000 C) are ~ 10^ 
times stronger than the average magnetic held (i?av ~ 2 C). It is possible that even if the solar 
average photospheric magnetic held had been much lower, the magnetic held pressure inside a 
sunspot could still have reached the photospheric pressure, i.e., the ratio would have been much 
larger than ~ 10^. Since convection seems to have a major role in concentrating the magnetic held 
to form sunspots (Priest 1987, chap. 8), and ACB stars have strong convection, I assume the ratio 
of spots to average magnetic held to be larger than that in the sun by an order of magnitude. This 
assumption will have to be tested when more detailed dynamo models become available. The 
maximum value of magnetic pressure that can be attained inside stellar spots is therefore 


T^spot — 


SCVl 


(7) 


with f] ~ 10*^. 

(2) For a spot to be signihcantly cooler than its surroundings I assume that the magnetic pressure 
is of the order of the photospheric pressure 


r2 

-^spot 

Stt 


~ R 
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( 8 ) 


(3) The pressure of the photosphere is taken to be (e.g., Kippenhahn & Weigert 1990, §10.2) 

_ 2 GM 1 


( 9 ) 
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where k ~ 3 x 10 ^ cm^ g ^ is the opacity at the surface of the AGB star. This matches the 
photospheric pressure of numerical models quite well (e.g., Soker 1992). 

Using the three equations (7)-(9), which serve as the basic assumptions of the model, gives for 
the average magnetic intensity required to form AGB stellar spots 


-Bav >4 X 10 


-3 (JlV'"( 


\lMs 


R* 


V300B(7 
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( 10 ) 


where = k/3 x 10“^ cm^ g“^. Taking Bav from equation (10) in the expressions for the magnetic 
energy luminosity Eb = 47ri?^u^Bav/87r, and using the kinetic energy luminosity Ek = Mwv‘^/2, 
give the magnetic activity required for the formation of AGB stellar spots 
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(11) 


As mentioned earlier, there is no basic theory to predict the magnetic activity of stars. Soker 
& Harpaz (1992) estimated, through a phenomenological study, the magnetic activity of rotating 
AGB stars. They assume that the amplification time of the magnetic field is equal to the rotation 
period. Using their equations (2.7) and (2.8), with the assumption of amplification time equal to 
the rotational period, with minor modifications of scalings, the ratio of magnetic to wind energy 
loss is 


^ ^ 4/ ^ Y f ^ ^ / M \ / Menv \ / Mw \ ^ / Vw 

\Ek)^^^^^~ ylO-^WKepy V3OOM0y \IMq ) \0.1Mq J Mq yv-^ ) UO km s”! 

( 12 ) 

By comparing equation (11) with equation (12), it turns that according to the dynamo activity 
estimated by Soker & Harpaz (1992), the required AGB stellar angular velocity is ~ 10“^ times the 
Keplerian surface velocity. The dynamo activity derived by Soker & Harpaz (1992), however, follows 
ideas developed for the solar dynamo, where the rotation velocity is faster than the convective 
velocity, opposite to rotating AGB stars. It is therefore possible that the amplification in the 
strong convective envelope of AGB stars is more efficient. I propose that the model presented here 
is effective for 
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(13) 


4. DISCUSSION 


Several points emerge from the analysis presented in the previous section. 

(1) By comparing equation (13), for the required angular velocity, with equation (2), for the angular 
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velocity of singly evolved AGB stars, it turns that single stars may possess high enough angular 
momentum to account for axisymmetrical mass loss, according to the proposed model. However, 
this requires that the surface angular velocity on the main sequence be ^ 5% of the Keplerian 
angular velocity. The sun, for example, is on the border; it rotates very slowly, but the fraction of 
the envelope that will be retained on the upper AGB, when the envelope mass is, say, 0.05Mq, is 
larger than the scaling used in equation (2). Other models, such as those described in §2, require 
substantial spin-up by a companion. This led me (Soker 1996; 1997; and references therein) to 
propose that planets are commonly present around stars. However, if in the next few years the 
results of the intensive planet search projects is that -C 50% of all stars have planets, then there 
will be a need for a model of efficient axisymmetrical mass loss for singly evolved stars. In addition, 
the mechanism should account for the increase in the degree of asymmetry toward the termination 
of the AGB evolution (see §1). The mechanism proposed in the present study may be such a 
mechanism, as is the mechanism of mode-switch to nonradial oscillations, proposed by Soker & 
Harpaz (1992). It is hard to tell from this preliminary study if the mechanism proposed here can 
be applied to a large fraction of singly evolved stars, since most main sequence stars may not rotate 
fast enough. However, the mechanism is quite effective if a planet companion of mass ^ O.lMjupiter 
enters the envelope at late stages of the AGB or the RGB and spins-up the envelope. 

2) From equation (2) we see that the angular velocity decreases rapidly as the envelope mass 
decreases toward the termination of the AGB. In the present model this decrease is more than 
compensated for by the increase of the vulnerability of dust formation and photospheric conditions 
to the magnetic activity (point 6 below). 

3) From equation (6) and the requirement a ^ 10“^, we find for the magnetic model of 
Ghevalier & Luo (1994) and Garcia-Segura (1997), Es/Ek ^ (^;rot/0.01^;^^,)“^. It should be noted 
that equation (11) refers to AGB stars, while the model of GL refers to post-AGB stars, which are 
much hotter and smaller, but are left with a small envelope mass, and are expected to rotate very 
slowly according to equation (2). It is clear from comparing equation (11) with the value given 
above that the presently proposed model is much less demanding on the magnetic activity and/or 
rotation velocity on the AGB and beyond, compared with the models described in §2. 

4) The required magnetic activity Eb (eq. 10) does not depend on the mass loss rate. 

5) Because the magnetic energy released through the photosphere is much below both the 
kinetic and thermal energy carried by the wind, the magnetic activity will not heat the region 
above the photosphere, except perhaps in localized regions where the magnetic energy becomes 
extremely strong. This is the opposite situation to that in the solar corona, where the magnetic 
energy is larger than the wind’s kinetic energy, Eb/Ei^ ~ 500. 

6) How is the increase in the asymmetry degree with the increase of the mass loss rate toward 
the end of the AGB explained in the present model? The answer suffers from the lack of a theory 
to predict the formation of cool spots and the formation of dust above these cool spots. But let me 
try to propose two plausible explanations. 
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6.1) As the star ascends the AGB it expands and loses mass. Both processes cause a rapid decrease 
of the density and pressure in the envelope below the photosphere (e.g., Soker 1992). This is the 
region where the magnetic field behavior determines the appearance of the magnetic field in the 
photosphere. Therefore, the low density and pressure may result in a much more violent magnetic 
behavior, leading to the appearance of larger and more frequent stellar cool spots. A very low 
envelope density and pressure exist in extended envelopes formed by late final helium shell flash, 
such as proposed for R Coronae Borealis stars and the PN A 30 (§3.1). In these types of objects 
dust formation may be facilitated by the high metalicity as well. 

6.2) For the mechanism proposed here to be effective, dust formation should occur close to the spot, 
since the temperature difference between the spot and its surroundings is smoothed and vanishes as 
the distance from the stellar surface increases (Frank 1995). As the star gets cooler, as is the case 
for stars ascending the AGB, the dust forms closer to the star, and hence the cool spots influence 
dust formation more strongly. The transition to a carbon-rich AGB star on the upper AGB may 
also make dust formation easier (Frank 1995). 

7) The higher mass loss rate during the superwind phase is explained both by the increase 
of the density scale height (Bedijn 1988; Bowen & Wilson 1991; §2.1), which influences the global 
mass loss rate, and by the increase of the magnetic activity, with increasing activity toward the 
equatorial plane. Therefore, both the axisymmetrical mass loss, which results from the magnetic 
activity, and the total mass loss rate increase as the star evolves along the AGB. The existence of 
spherical PNs which have no superwind, ~ 75% of all spherical PNs, has two plausible explanations 
in the frame of the proposed model. The first is that the progenitors of these PNs lost their entire 
envelope during the early AGB. The second is that the magnetic activity is responsible for the 
superwind in some cases, and the progenitors of the spherical PNs were rotating too slowly to 
support magnetic activity. 

8) The solar magnetic activity has a cycle of 11 years period. What will be the consequences, in 
the proposed model, if rotating AGB stars have magnetic activity cycle? It is impossible to predict 
the period of such a cycle because of the lack of a detailed dynamo model. But if exist, it will, 
according to the proposed model, cause oscillations in the mass loss rate. I would like to speculate 
that the arcs, or shells, found recently in several PNs (e.g., GRL 2688 [Egg Nebula], Sahai et al. 
1998a,b; IRAS 17150-3224, Kwok, Su, & Hrivnak 1998) may result from such magnetic activity 
cycle. In IRAS 17150-3224 the arcs are concentric, nearly circular (Kwok et al. 1998), while in 
GRL 2688 the arcs show departure from circularity, and a large fraction of them span only small 
angles (Sahai et al. 1998b). They arcs represent almost periodic enhancement in mass loss rate, 
(Kwok et al. 1998), by a factor of at least 2 (Sahai et al. 1998a). The periods are a few hundred 
years. A magnetic activity cycle, similar to the solar cycle, but of a few hundred years period, 
may explain the almost periodic, and somewhat sporadic in direction, mass loss that formed the 
shells. A possible problem to the proposed model is that the shells are almost spherical, though 
some differences between the equatorial and polar direction exist. Other models for the formation 
of the arcs (e.g., Harpaz, Rappaport & Soker 1997) are discussed by Sahai et al. (1998b). 
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5. SUMMARY 


The main goal of the paper was to propose a mechanism for axisymmetrical mass loss on the 
AGB that: (a) increases as the star evolves along the AGB, and (b) operates for slowly rotating 
AGB stars (having angular velocity in the range of lO^^WKep ^ ^ lO^^WKep, where WKep is 

the equatorial Keplerian angular velocity). Such angular velocities could be gained from a planet 
companion of mass ^ O.lMjupiter, which deposits its orbital angular momentum to the envelope 
during the AGB phase or even much earlier during the RGB phase, or even from single stars which 
are fast rotators on the main sequence. The model is built to explain elliptical PNs, but not the 
more extremely asymmetrical bipolar PNs, which are thought to be formed from stellar binary 
systems. 

The proposed model assumes that dynamo magnetic activity results in the formation of cool 
spots, above which dust forms much easily. The enhanced magnetic activity toward the equator 
results in a higher dust formation rate there, hence higher mass loss rate. As the star ascends the 
AGB, both the mass loss rate and magnetic activity increase rapidly. 

Future observations should search for magnetic fields around AGB stars, e.g., in masers, and 
for sporadic mass loss episodes, such as in R Goronae Borealis stars. An extremely important issue, 
which is relevant to all models of axisymmetrical mass loss, is the detection of rotations in AGB 
stars. On the theoretical side, it will be necessary to extend the work of Frank (1995), i.e., calculate 
the density and temperature profile above AGB stellar spots, and from that the dust formation 
rate. Gool spots formed by large convective elements (Schwarzschild 1975) will also have to be 
considered. These are difficult tasks because of the strong pulsations in AGB stars. An even more 
difficult task is to build a dynamo theory for slowly rotating AGB stars. 
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